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Abstract 

A simple explanation of the Karman vortex street is suggested on 
the basis of the two-temperature canonical distribution for inviscid 
two-dimensional flows in Navier-Stokes fluids or guiding-center plasmas. 
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The Kill man vortex street develops ns a final state for a variety 

of unstable shear flow profiles* In two-d linens ional Navier-Ztoken 
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fluids or electrostatic guiding center plasmas. The vortex street 

is an old and familiar experimental phenomenon^ which has apparently 

never been explained theoretically. Though painstaking and accurate 

numerical solutions of the Navicr-Stokes equations have followed the 

development of initially unstable shear flow profiles into vortex streets 

(see particularly Zabusky and DcemS , there is apparently no simply 

accessible theoretical picture of why the two rows of oppositely-signed 

vortices arrange themselves as they do in the characteristically staggered 

arrangement across the street. The purpose of this letter is to suggest 

a simple explanation of the vortex st.eet on a statistical-mechanical 

basis. 

Our slatting point is a generalised application of the two-temperature 
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Ci inical distribution of K-iichnan * for inviscld, two-dimensional, 

Navier-Stokes flows. The coordinates of the phase space over which this 

ensemble is defined are the real and imaginary parts of the velocity 
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field Fourier coefficients in a truncated Fourier representation. The 
canonical distribution varies as cxp(-aE - Bi2) , where E and ft are the 
mean energy and "enstrophy" (mean square vorticity) densities, a * and 
B * are two temperatures, one for energy and one for enstrophy. The 
relaxation of arbitrary initial conditions to final states which have 
properties that are accurately predicted by this ensemble has been given 
extensive numerical investigation by Seyler et nl . *** Relaxation to this 
equilibrium state now appears to be beyond dispute. 
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One of the more interesting predictions of the two-tcmpctature 

canonical distribution is for the expectation value of the kinetic 
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energy density associated with the k th Fourier component: 

<|u(k)| 2 > - (a + PkV 1 . (1) 


The angular brackets < > denote an average over the two-temperature 
canonical distribution. The reciprocal temperatures a and 8 are in 
effect Lagrange multipliers whose values are determined by the requirements 
that Eq. (1) lead to the correct initial values of E and ft: 


E - L\ (a + $k 2 ) -1 
*k 

( 2 ) 

ft - k 2 (a + 3k 2 ) ~ 1 
k 


The interesting thing about Eq. (1) is that qualitatively different 
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spectra and flow patterns ’ are associated with different regimes of a 

and 3 . Either a or 3 (but not both) can be negative, depending upon the 

ratio ft/ E. In particular, for either large k or low values of ft/E, 

max 

2 

a<0, 3>0, and a + 3k is only slightly greater than zero at the minimum 
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allowed values of k . In these cases, the gross features of the macro- 
scopic turbulent field are dominated by the first few Fourier coefficients 
(longest wavelength modes) and become relatively insensitive to the higher 
wave-number contributions. The flow consists of a large pair of counter- 
rotating vortices that fill up the box inside which the Fourier analysis 
is performed. ^ This basic pattern had previously been noted in numerical 
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solutions for the finite viscosity case by Deem and Zakusky and 
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Tnppert and Hardin, and for the discrete linc-vortox representation 
by Joyce and Montgomery. ^ 

The point to be made now is that the two-temperature canonical 

distribution has a wider applicability than simply to the case of periodic 

boundary conditions In a square box. For more elaborate geometrical 

arrangements in which fully developed two-dimensional turbulence is 

expected, another set of eigenfunctions than simply all the plane waves, 

exp(ik • x) , may be appropriate. Statistical-mechanical analyses of hydro- 

dynamic situations of a Tre general character using an arbitrary set of 
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basis functions have bei suggested by Thompson. 

We discuss a situation qualitatively similar to one considered by 
Zabusky and Deem,^ though detailed comparisons are not in order (we omit 
viscosity). An unstable shear flow profile with velocity parallel to 
the x-axis is assumed to develop into a vortex street between y = -L /2 
and y = +1,^/2. The fluid velocity u(x) has only x and y components and 
is independent of the z coordinate. The vorticity p(x) = Vxu = p(x,y)e 

Z 
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lies normal to the xy plane. The stream function ij)(x,y) obeys V <j> = -p 
and is related to the fluid velocity by u = V<f) x e . (Each realization 
of the ensemble will of course be time-dependent, but we suppress the 
time arguments.) We assume periodicity in x with a period L . We require 


that there be rigid frictionless walls at y = +L /2, so that e • ( V d> x e ) 

— y y z 

vanishes there. The effect of the boundaries can be replaced, using the 
method of images, by an infinite periodic array of image vorticity with 
period 21.^ in the y-direction. The vorticity elements outside the basic 
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box of width L arc chosen so that equal ana opposite elements are always 
connected by straight lines which are perpendicularly bisected by the 
walls. A situation studied by Zubusky and Deem was for initial velocity 
u “ u° e , where u° » U - U exp(-y^/A^) , with L »2A, so that the 
shear flow velocity was essentially zero at the walls. It is convenient 
to work in a coordinate system in which the final vortex street is at 
rest, which implies a non-zero overage flow velocity in the x-dlrection. 
Since the absolute square of each Fourier coefficient except k = 0 is 
invariant to a Galilean transformation anyway, the mean square vorticity 
distribution can easily be transformed from one coordinate system to 
another. The k “ 0 component of the fluid velocity is determined strictly 
by momentum conservation, and is not involved in a discussion of the 


interactions of the k ^ 0 modes. 

The most general representation of the stream function under the above 
periodicities is (again suppressing the time arguments and omitting the 
part associated with tho uniform translation): 


<Kx,y) » /L tA x <V n v ) sin k x x 

n ,n 3 

x’ y 


+ B (n ,n ) cos k x] [A (n ,n ) sin k y 
x x’ y x 1 y x* y y J 


+ B (n ,n ) cos k y] (3) 

y x* y y-' ' 

where k = 2ir n /h , n =0, 1, 2, 3, ..., and k “ 2u n /2 L , 

n^ = 0, 1, 2, 3, ... . However, the geometrical requirements of the 

situation demand that certain of the A and B coefficients vanish. First, 

(n ,n ) = (0,0) does not become involved in the dynamics and may be omitted, 
x y 
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More Importantly, the condition that e • (e x V<jO * 0 at y ■ + L /2 

y z y 

demands that for n 4 0, A (n ,n ) » 0 for n odd or zero, and that for 

x y x y ■ y > 

n 4 0, U (n ,n ) » 0 for n even or zero. For n » 0, the presence of 
x y x y y x * 

the walls Imposes no additional condition, but the physical requirement 

that the average x velocity at the upper and lower boundaries of the box 

be equal requires that By(0,l) “ 0. 

For the present situation, the prediction of the r.,can modal energy 

given by Eq. (1) is still correct in the frame with zero total momentum, 

as long os one keeps in mind that the space of eigenfunctions has been 

somewhat pruned. For large values of E/f2 (and since in the Zabusky-Deem 
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example, this ratio is just A as long as A«Ly, this can be as large 

2 

as we like), we are in the regime a<0, £S>0, and a + 3 k ^ VGr ^ small and 

►V . 

positive, so that only the lowest k-modes dominate. Since k jt 0 for the 

modes described by Eq. (1), the expectation values are the same in any 

frame, including the one in which the contours of constant vorticity are 

expected not to translate. The problem reduces itself o that of 

enumerating the allowed modes with non-vanishing Fourier coefficients 

2 2 

and the lowest values of k + k . Since the (n ,n ) - (1,0) is ruled 

x y x y 

out, the first two modes clearly are (n ,n ) - (0,1) and (1,1). Since 

x y 

B (0,1) “ 0, only the sinTiy/L term is permitted. To .’no extent that 
2 2 

43ky»a + 3k y >0, values of n^>l need not be considered. The energy 

associated with the (1,1) mode is comparable to that of (0,1) to the 

2 2 2 2 
extent that 3^ x is of the order of a + 3 k^ (this implies k»k^, or 

L »L ). The energy of the (1,2) mode is not totally negligible compared 
x y 

to that of the (0,1) mode under these circumstances (a factor of about 

2 2 

1/6 follows from assuming a + 3 ^ » 3 k ) , but for purposes of determining 

y x 
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the qualitative behavior, signif leant Insight can be gained by estimating 

the flow field solely on the basis of the contributions from the (0,1) and 

» 

(1,1) modes plus the uniform translation associated with the (0,0) mode. 
Znbusky and Deem found in fact that the significant fraction of the energy 
was contained in their first three modes . ^ 

Irrelevant phase factors, not provided by the canonical distribution, 
are associated with each realization of the ensemble and must be assumed 
in order to draw typical Clow fields. For illustrative purposes, there 
is no reason to assume they are anything but zero. For a first approxi- 
mation, then, we may assume for <p the approximate form (omitting the 
uniform translation) : 

<J)(x,y) = a sin (iry/L y ) 

+ b cos (iry/Ly) cos (2tt x/L^) , (4) 

where b/a is a number comparable to, but somewhat less than, unity. 

Figure 1 is a schematic representation of the way the maxima and minima 
of the simple approximate form (4) for <J) combine and interfere to produce 
the vortex street configuration with periodicity length L in the 

A 

x-directiun. Addition of higher terms of the series will modify the 
details of the contours of constant <J>, but cannot change the gross confi- 
guration. '/.t is therefore suggested that the two-temperature canonical 
distribution explains in an exceedingly simple way the basic structure of 


the Kdrmdn vortex street. 
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The intention is to carry out future numerical Invisdd solutions 

to the 2D huvicr-Stokcs equations in the present geometry, watch the 

♦ 

vortex street develop, and attempt a detailed comparison of the time- 
averaged modal energies with the spectral predictions of Eq. (1). 

The author thanks Drs. P. D. Thompson and C. E. Leith each for a 
helpful conversation. Part of this work was performed at the University 
of Iowa under NASA Grant NGL-16-001-043. 
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FIGURE CAPTION 


Fig* 1 Schematic arrangement of the maxima and minima of the stream 

function (the uniform translation omitted) nnd/or the vorticity. 
The upper panel in for the first term in Eq. (4), the middle 
panel for the second, and the bottom panel for the sum. Without 
loss of generality we may assume a>b>0. The values chosen, for 
illustrative purposes only, arc a ® 1.0, b “ 0.8, 






